Universal and non-universal amplitude ratios for scaling corrections on Ising strips 
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We consider strips of Ising spins at criticality. For strips of width A'^ sites, subdominant (additive) 
finite-size corrections to scaling are assumed to be of the form a^/N^ for the free energy, and h^/N^ 
for inverse correlation length, with integer values of k. We investigate the set {0^,6^} [k > 2) 
by exact evaluation and numerical transfer-matrix diagonalization techniques, and their changes 
upon varying anisotropy of couplings, spin quantum number S, and (finite) interaction range, in 
all cases for both periodic (PBC) and free (FBC) boundary conditions across the strip. We find 
that the coefficient ratios bk/a^ remain constant upon varying coupling anisotropy for S = 1/2 and 
first-neighbor couplings, for both PBC and FBC (albeit at distinct values in either case). Such 
apparently universal behavior is not maintained upon changes in S or interaction range. 

PACS numbers: 64.60.De,75.10.Hk,05.70.Fli 



I. INTRODUCTION 

In this paper we investigate corrections to scaling in 
critical Ising systems on a strip geometry. Consider a 
square lattice with N lines and M columns, in the limit 
M ^ 00. Other two-dimensional lattices, such as trian- 
gular or honeycomb, can be brought into a square-like 
shape, by suitable bond additions or deletions. From the 
largest (Ao) and second- largest (Ai) eigenvalues of the 
column-to-column transfer- matrix (TM) , one obtains the 
free energy per spin, /jv (in units oi ksT), and spin-spin 
correlation length ^jv, via 



iV/jv=ClnAo 



N 



/■ 1 ^0 



(1) 



The factor C, is unity for the square lattice and, in tri- 
angular and honeycomb geometries (also for the square 
lattice when the TM progresses along the diagonal [2,l3l)> 
corrects for the fact that the physical length added upon 
each application of the TM differs from one lattice spac- 
ing In all cases of interest here, i.e., ferromagnetic 
systems, Aq and Ai are both real and positive. 

At the critical point Tc where a second-order transition 
takes place, conformal invariance Q gives the following 
relations regarding universal quantities c, the conformal 
anomaly (6|, and the spin scaling dimension xi 0|: 
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In Eq. ([2]), where c = 1/2 for models in the Ising univer- 
sality class, /surf = for strips with periodic boundary 
conditions (PBC) across, and non-zero for free (FBC) 



or fixed BCs; a = ^ for PBC, and ^ for FBC In 
Eq. ([3]), where the exponent xi for the Ising universality 
class is = i in the bulk, and xf — ^ for the ordinary 
surface transition, one has xi ^ x\, (3 ^ 2 ior PBC, and 
xi = xf, /3 = 1 for FBC [7]. 

Since Eqs. ([2|) and ([3]) are expected to be exact only 
asymptotically, it is of interest to develop a systematic 
understanding of the corresponding finite- iV corrections. 
We write: 



A^(/jV-/oo)-2/3urf =^ 



k=l 



ah 
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where oi — aire, hi — I3t:xi. Assuming only integer pow- 
ers of N^^ in Eqs. ^ and ([5]) is believed to be warranted 
as long as one is dealing with models in the Ising univer- 
sality class [H . We revisit this assumption in Sections IIVI 
and |V] below. Our task here will be to learn as much as 
possible about the coefficients {ak, bk}, k > 2, as well as 
(for reasons explained below) their ratios bk/ak- We are 
interested in their respective universality, or lack thereof, 
upon changes in boundary conditions, degree of spatial 
anisotropy of interactions, spin quantum number S, and 
(finite) interaction range. We restrict ourselves to the 
square lattice. 

In Section |TT] we investigate 5* = 1/2 strips with PBC, 
first-neighbor interactions, and varying degrees of spatial 
anisotropy; in Sec. Illli we examine systems with FBC, 
again with varying anisotropy; Sec. IIVI deals with the 
spin-1 case, and isotropic couplings only; in Sec. |V] we 
return to ~ 1/2 and introduce next-nearest-neighbor 
couplings (keeping to isotropic interactions). Finally, in 
Sec. I VI I concluding remarks are made. 
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II. PERIODIC BOUNDARY CONDITIONS 

A. Preliminaries; isotropic systems 

We recall results for strips cut along the x direction, 
with N lines and M — >■ c» columns, and PBC across. All 
eigenvalues of the TM can be written in closed form 
With Ki = Ji/ksT being the interactions respectively 
along X {i — 1) and y {i — 2), Aq and Ai are: 



InAo - iA^ln(2sinh2i^i) = - ^ 

1 1 

InAi A^ln(2sinh2ii'i) = - ^ 72^ 



l2r+l , 



(6) 



(7) 



r=0 



where 



cosh7r = cosh2i^^ cosh2i^2~sinh2iirj' sinh2iir2 cosw^ ; 

(8) 

the dual couplings K* are defined by tanhX* = 
exp{—2Ki), and the allowed frequencies are ojr = rn/N. 

With Si = smh2Ki, one has si S2 = 1 at the critical 
temperature where the system is self-dual, and Eq. ^ 
becomes: 



1 



cosh7r = 1 H — 2 ~ coswr) (T = Tc) 



(9) 



For isotropic systems, si = S2 = 1 at criticality. In 
this case, the sums in Eqs. © and ([7]) were tackled Q 
by applying the extended Euler-MacLaurin summation 
formula fm,^- 



N-l ^ „b 

^ F{a + nh + ah) = - / F{x) dx + 



fc=i 

where h = {b - a)/N, F'-^^x) is the j-th derivative 
of < a < 1, and the Bk{a) are the periodic 

Bernoulli polynomials (related to the BernoulUi numbers, 
denoted simply by Bk, by Bk = Bk{Q)). 

It was found that only odd powers of N~-^, i.e., k — 
2j — 1, j > 1 occur in Eqs. (|4]) and (O; this can be 
traced back to the fact that the Bernoulli numbers _B,„ 
obey i?2m-i = (m > 1). Also, relatively simple closed- 
form expressions were derived for all and bk ■ Such ex- 
pressions reproduce previously-known exact results (for 
ai 61 Q, and 63 [13 )j and are in very good agree- 
ment with numerically-obtained ones (l3l |. Furthermore, 
although the coefficients themselves are non-universal 
(upon changing lattice structure, or considering quan- 
tum Ising chains [l3 - [l6| instead of their two-dimensional 
classical counterparts) , their ratio is found to remain con- 
stant upon the same set of changes [9]: 



bk_ 

a-k 



2fe+i 



2*^ - 1 



(fc = 2j-l , j> 1) [PBC]. (11) 



It should be noted that when one considers the TM 
running along the diagonal of the square lattice (as in 
Refs. [2 andy), one gets for isotropic systems with PBC 
the same value for the ratio bk/ak as in Eq. pTjl . Further- 
more, the coefficients themselves have the same absolute 
value as those found with the TM along x; only, they 
alternate in sign: a^, bk < for /c = 1, positive for k ^ 2 
etc [13. 



B. Anisotropic systems with PBC 

With K1/K2 = R ^ 1, one gets fis', ^ the corre- 
sponding forms of Eqs. ^ and Q [specializing to Ising 
spins on strips with PBC ] as: 



lim (/a,(,) - /oo) 



1 TT 

V,, 12 



(12) 
(13) 



where /Ar(i) and are, respectively, free energy and 

correlation length at criticality. both calculated by it- 
eratin g th e TM along the direction with couplings Ki. 
Note [13I that /oo in Eq. also depends on R . 

si = sinh2ii'c is the solution of sinh2iirc smh2RKc = 1. 
As noted in Ref. 1^ Eq. Q can be rewritten as: 



7(w) = 2 ln(u + Vl + m2) , 



1 . LJ 

— sm — 

si 2 



(14) 



In this form, it is immediate to see that anisotropy brings 
about a simple rescaling of the argument in the sums of 
Eqs. (ini) and ([7]). Furthermore, in the Euler-Maclaurin 
formula, 7(w) only occurs through its derivatives of n-th 
order at the endpoints uj = and tt, which satisfy 
7("-)(7r) = -^(")(0), see Eq. This is enough to guar- 

antee that any coefhcient ak{R) [bk{R)] will differ from 
its isotropic counterpart afe(l) [6^(1)] by a multiplicative 
correction, 5fc(si). Thus, it is predicted in Ref. that the 
ratios given in Eq. (|lip will remain unchanged. In this 
context, Eqs. and reflect the (easily checkable) 
fact that 51 (si) = Si^, where for Eq. (fT5|) one also uses 
si S2 = 1 at criticality. 

In order to test the robustness of the theoretical frame- 
work just expounded, we evaluated the third-order cor- 
rection. This is done by replacing the argument of 
Eqs. (ini) and ([7]) by its generalized form, Eq. (fn| . and 
following the corresponding effects on the term in 
Eq. (fTU|) , which arise from the third-order derivatives in- 
dicated there. One finds: 



53(si) 



1 

2^'i 



2(S1 



Si ' 



2siti 



(15) 



where ti = tanh 27-^1. 

We numerically calculated /^r and from 
Eqs. ([1]), ([6]), and (O for assorted values of R, and 
N = 10/, / = 2, 3... 30. The resulting sequences were 
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Figure 1. (Color online) Points: results for a3{R), 63(ii) de- 
fined in Eqs. p6p and (|17|) . from fits of /at, data evalu- 
ated for strips of widths 100 < N < 300 and PBC across, via 
Eqs. (O, and ([T]). Uncertainties are smaller than symbol 
sizes. Full lines show (73(51) defined in Eq. (llSp . multiplied 
respectively by a3(l) [blue] and 63(1) [red] (see text). 



It is interesting to consider the above results for 
R 00. It is known [H, [H, [2l|, [Hi that the [zero- 
temperaturel quantum Ising chain (QIC) in a trans- 
verse field [Mj has a correspondence with this extreme 
anisotropic limit, via: /oo £-01 Ei — Eq etc, 

where the Ei are the energy levels of the quantum sys- 
tem. In Ref. Q the energy spectrum of the QIC with PBC 
was studied directly with help of the Euler-Maclaurin 
formula, and the corresponding ratio bk/ak was found 
to obey Eq. ((TT|). The latter result can also be ex- 
tracted from the exact expressions Eqs. (17a) and (18a) 
of Ref. 21. 

While the limit given in Eq. is preserved as i? — > 
cxD, the exponential divergence of the higher-order terms 
is not cancelled: 



1 ' = h - (si H ) 

AN 4 ^ si 



(18) 



A similar effect (with the factor 1/si) is already obvi- 
ous in Eq. (|12[) . In summary, although coefficient ratios 
bk/cLk are preserved as i? — > cxd, each term of the Euler- 
Maclaurin expansion for the two-dimensional Ising model 
with i? ^ 1 is translated into its counterpart of the cor- 
responding expansion for the QIC by means of a distinct 
anisotropy factor. 



III. FREE BOUNDARY CONDITIONS 



adjusted to: 



In{R) = /oo(i?) + 



1 



03 (i?) , a^iR) 



si 127V2 
bm 



iV4 

b5iR) 
N5 



Are 



(16) 
(17) 



where /oo, {ofe}, and {bk} are adjustable parameters. It 
is important to keep the next-higher-order terms as and 
65 in the truncated expansions above, in order to im- 
prove stability for the quantities a^, and 63 which are the 
main focus of interest here. The optimum range of N, 
large enough for higher-order terms to have negligible 
influence, but not so large as to compromise the numer- 
ical accuracy of fits (since this depends crucially on dif- 
ferences between finite- iV estimates of and ■Civ^) 
found to be 100 < iV < 300. In Fig. [1] we show 03 (i?) and 
63 (i?), fitted via Eqs. and ^T7\ . for several values of R 
spanning four orders of magnitude. The continuous lines 
depict Eq. ([T5|) . multiplied respectively by the isotropic 
values 03(1) = 77r3/1440, 63(1) = 7:^/96 SEIll. The 
agreement is perfect, except for 03 (i?) at i? > 30 where 
reasonable convergence was only obtained upon adding 
the next higher-order term, a'j{R) /N^ , in Eq. (fTO)) . 

Our results provide direct numerical evidence that the 
Euler-Maclaurin scheme used in Ref. is indeed appli- 
cable to Ising systems with PBC and any finite degree 
of (ferromagnetic) anisotropy. Similar conclusions were 
drawn for the anisotropic Ising model with Braskamp- 
Kunz boundary conditions [20,] . 



A. Isotropic systems 

The eigenvalue spectrum of the TM has been ob- 
tained |23l - [25j for Ising S — 1/2 strips with nearest- 
neighbor couplings and free boundary conditions (FBC) 
across. For a strip of width N sites, one has: 



1 ^ 

In A„, = 2 X! ±7(wi) , m = 0, ■ 
1=1 



(19) 



where the ± combinations run through all 2^ possibili- 
ties. A regular background term, iA^In(2sinh2iiri) [see 
Eqs. (O and © ], has been omitted. With ah the 7(0;;) 
real and positive for this case [23j . 



N 



InAo 



N 



n=l 



UJ„ 



(20) 
(21) 



where uii corresponds to the smallest 7. The relationship 
between the 7 and the allowed frequencies uii is given 
by 0: 

cosh 7 — cosh 2_fi'i cosh 2K2 — sinh 2Ki sinh 2K2 cos lu ; 

(22) 



sinh 7 cos S* = sinh 2K2 cosh 2_fi'i — 
— cosh2_ft'2 smh2Kl cos a; ; 



(23) 
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Figure 2. (Color online) Ising strips with FBC and isotropic 
interactions, at Tc'. differences tj„ — ijJn-i {n > 2) between 
consecutive solutions of Eqs. p6|) and (|27p . for strip of width 
N = 30 sites. 



Figure 3. (Color online) Ising strips with FBC and isotropic 
interactions, at Tc'. illustrating the convergence of N{ujn —(^55) 
towards f{n/N) with increasing A'^; see Eq. (|28|) . The inset is 
a blowup of the central section of the main figure. 



sina;/sinh7 = sin(5*/sinh2Xi . (24) 

From Eqs. (EH), (gSl), and ([M]), one gets at the critical 
point, where si S2 — 1: 

cosh 7 = 1 + ^(1 — cosoj) ; (25) 

•n 

tan<5*=ti-^^^, (26) 
(1 ~ cos w) 

again with ti = tanh2JCi. From Eq. (P5|) . the smallest 
7 corresponds to the lowest allowed cj. Note also that 
Eq. ((25)) is identical in form to Eq. ([9]), so it can also be 
rewritten as Eq. Finally, the allowed frequencies 

Lu„ can be determined from Eq. (|26p combined with the 
quantization condition [23|: 

^^N^ ^ ^^^S' ^ < W < TT , (27) 

by eliminating the auxiliary angle S* . 

For the remainder of this Subsection, we shall consider 
only isotropic systems (A'l = K2), thus si = 1, ti = 
l/\/2 in Eqs. 1^ and 

The resulting frequencies are not equally spaced, as il- 
lustrated in Fig. [21 So the Euler-Maclaurin formula can- 
not be used in the same way as in Ref. !). to calculate 
the free energy from Eq. (EUl) . However, we show in the 
following that one can still make adaptations and extract 
some useful information. We found that for large N the 
LOn approach the form: 

(28) 



As shown in Fig.|3l f{u) is a smoothly varying function 
of u = n/N. One has: /(O) = 0, lim„_,i f{u) = -tt/2. 
Both limits can be understood by examination of the 
graphical solutions of Eqs. and (gTJ) [Hi. The resid- 
ual iV-dependence of N{ujn — ^n) highlighted in the 
inset of Fig. [31 This can be accounted for by an addi- 
tive correction of the form (l/N) f2{n/N); /2 is nearly 
constant, varying smoothly between sa 1.1 and 1.3 for 
< n/N < 1. We thus write: 

00 

^{u^) = j{coO) + J2 — gO){u:'J, (29) 

with 

gW(c.°)^{/(c.)7(^H^)}^^ , 

ff(2)(^°)^|/2(w)7WM + ^/'M7^'H^)} ^ ■••(30) 

where 7*^™^ = d"^^/doj"^, and the arguments of / and 
/2 have been straightforwardly changed. Eq. ([^0)1 then 
becomes: 

N N 00 

lnAo = i^7(-°) + ^E5^'^(-«) + ----E^- 

71—1 n—1 j—0 

(31) 

So, each term (of order j > 0, with g^^\oj^) = 7(0;°)) 
of the Taylor expansion indicated in Eq. gives rise 
to a sum Sj of N terms, each of the latter evaluated at 
uj = Lu^ {1 < n < N), i.e., at equally spaced intervals. 
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We investigated the feasibihty of applying the Euler- 
Maclaurin formula, Eq. (fTU| . to each Sj, with x = to, 
h — tt/N, a — 0, b ^ tt, a — 1/2, so that the result would 
be of the form Sj = E^-i a- /iV(*+^) [where i = -1 
corresponds to the integral in Eq. ] . Sj would then 
give contributions to InAo at all orders A^~(J+') i> —I. 
Note that = because Bi(l/2)=0 [13, \lM- How- 
ever, one would have to assume that the infinite sum 
implicit in each Taylor series commutes with the infinite 
sum present in each separate Euler-Maclaurin expansion 
(the form given in Eg. pop assumes that the remainder 
term vanishes; see, e.g., Ref . ITTh . Having in mind that 
the expansion parameter of the Taylor series and the sam- 
pling interval of the Euler-Maclaurin formula can be of 
the same order (tt/TV), it is doubtful that such commu- 
tation can be guaranteed. 

With these words of caution in mind, here we evaluate 
only a few of the lowest-order terms which would occur 
in such a calculational framework. 

We applied the Euler-Maclaurin formula to Sq in 
Eq. (|3ip. This differs from the sum indicated in its PBC 
counterpart, Eq. ([5]), in that the frequency spacing here 
is half that in the latter Equation. For the corresponding 
integral of Eq. pTI)l , this is compensated by the fact that 
the integration interval is cut in half as well, so from a'Li 
one reobtains the bulk result foo — (1/2) ln2 = (2G/7r), 
G = 0.915965594... (Catalan's constant) [J. For the 
terms of Eq. ([TU]) involving derivatives of the m-th order, 
the corresponding term in So has an extra factor 2~(^+'") 
relative to its PBC analogue One gets = 7r/48 as 
given by conformal invariance 0, aS = 0.00942 ... [9]. 

For 5*1, we evaluated / = g^^' {u) du using finite- 
N approximations for f{x) with 200 < < 2000, and 
extrapolating the resulting sequence against 1/A^. The 
final result is //27r = a^^ = -0.1817309(1), to be com- 
pared with 2/surf = -0.18173148 . . . In the compu- 
tation of higher-order terms, we ran into inconsistencies 
between results thus obtained, and those coming from di- 
rect numerical evaluation of the free energy via Eq. (|20p . 
We conjecture that these difficulties stem from the con- 
ceptual problems in interchanging the order of infinite 
sums, referred to above. 

As regards the correlation length, from Eq. (PT|) above, 
and combining Eqs. ^ and ([Ti)) . one has for the finite- 
A^ estimate rjf (N) of the decay-of-correlations exponent 
7^1=2x1: 



2N . 



TT 



y 



coswi . (32) 



By solving Eqs. ^ and (HTj) in the hmit uj ^ 0, 6* ^ 
tt/2, and consequently taking y 1 in Eq. (|32p . one gets: 



vliN) = 1 



1 1 



21 



TT 

48 



1 



+ OiN-^) . (33) 



According to Eq. both odd and even powers of A^^^ 

are predicted to arise in the expansion of for this 
case. For Ising systems with FBC, the occurrence of A^~^ 



corrections to finite- A' estimates of scaling powers was 
noted in Ref. [H. 

We evaluated Jn and £,n for A = lOj, j — 2, ... ,30, by 
numerically solving for the allowed frequencies and then 
plugging the results into Eq. and, finally, Eq. p^ . 

We fitted free-energy data for 100 < Af < 300 to a 
truncated form of Eq. with k < 4. After ensur- 
ing that known quantities were reproduced to good ac- 
curacy when allowed to vary freely, we fixed them at 
their known values, naraely foe = (1/2) In 2 -I- (2G/7r) flj; 
/3^,f = -0.0908657 ... [lai; ai = 7r/48, with the results: 
02 = -0.04616(2), 03 = 0.024(1), = 0.69(6). Note 
that 03 as given here differs from aj] evaluated from Sq 
above, in connection with Eq. pip. This is because 03 
gets additional contributions from higher-order sums Sk , 
k > (not calculated there). 

A fit of a subset (100 < A < 300) of the ryf (A) thus ob- 
tained to the form rif{N) = i]f + J2k=i ^'k A^*^ gave rjl — 
1 (±1 X 10^1°), b[ = -0.70710680(3), 6^ = 0.294388(7), 
6^ = 0.2274(7), b[ = -0.68(3). By keeping ryf, b[, b^ 
fixed at the respective values predicted in Eq. (j33p . we 
obtained b^ = 0.227972(6), b^ = -0.7013(6). The above 
results both confirm the predictions of Eq. for b[ 

and 62, and indicate that, in general, both even and odd 
powers of A^"-'^ occur in the expansion whose lowest-order 
terms are given in that Equation. 

We defer analysis of the ratios bk/ak thus obtained 
until the next Subsection, where anisotropic systems with 
FBC, and their connection to the QIC with free ends, are 
discussed. 



B. Anisotropic systems with FBC 

We first note that, even though Eq. is valid here, 
the arguments given immediately below it do not seem 
to cover the present case, since for FBC the w„ depend 
on anisotropy in the non-trivial way given in Eqs. (j26p 
and ([?f)) . Thus it is not obvious whether, e.g., Eq. 
still applies to the free energy here. 

We have directly examined the for varying 

anisotropics, and seen that their behavior is qualitatively 
similar to that for the isotropic case, depicted in Figs. [5] 
and[21 In particular, the limits /(O) — and /(I) = — 7r/2 
still hold [see the comments following Eq. (1^51) ]. 

By incorporating anisotropy into Eq. via Eq. P^ , 
one gets the generalized version of Eq. (|55)) : 



1 

si 



1 1 

2hN 



1 


'1 


7r2- 


^1 


2tf 


2 


48 


iV2 J 



(34) 

We numerically calculated Jm and from 
Eqs. dl]), (|^ . and (|^ for assorted values of R, and 
N = lOj, j = 2, 3 ... 30. Bearing in mind the FBC- 
adapted forms of Eqs. (HI]) and (HH) [ll[ll|, the resuhing 
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remain constant against any finite degree of (ferromag- 
netic) anisotropy. Note that fea/oa differs substantially 
from the PBC value 15/7 Q. It is remarkable that the 
free-energy coefficients depend on anisotropy in the same 
way as those for the correlation length. As stated in 
the first paragraph of this Section, this is not obviously 
granted at the outset. 

We consider the extreme anisotropic limit i? —> oo of 
Ising strips with FBC, and its connection to the QIC 
with FBC at both ends [ll,|2ll. In Ref.H^, the exact ex- 
pressions for ground-state energy and energy gaps of the 
QIC, given in Ref. [ij, were written as Euler-Maclaurin 
expansions; similarly to the PBC case, only odd powers 
of N^^ were found to occur in the corresponding forms 
of Eqs. dH) and (P. The counterpart to Eq. in this 
case was shown to be: 

- = M (fc^2j-l, j>l)[FBC]. (37) 

flfc (2'=-l)_Bfc+i 



Figure 4. (Color online) Points: results for {ak{R), bk{R)} 
defined in Eqs. (|35p and H36p [absolute values for k = 2], 
from fits of /jv, ^J^^ data evaluated for strips of widths 
100 < N < 300 and FBC across. Uncertainties are smaller 
than symbol sizes. Lines show anisotropy factors defined in 
Eq. (|34p . multiplied respectively by afe(l) [blue] and fefe(l) 
[red] [k = 2,3] (see text). 



sequences were adjusted to: 

NifN-U)^2U, + l-J^+±^J^, (35) 

^ fe=2 
fe=2 

where /oo, /surf (both of which, as well as /at, also depend 
on R), {flfe}, and {bk} are adjustable parameters. The 
b'f., defined in connection with Eqs. ([5^ and relate 
to the bk of Eq. ^ by bk = {n/2) b',,_^. As done in 
Section [III we keep the next-higher-order terms a4 and 
&4 in the truncated expansions above, in order to improve 
stability for the quantities 02, 03, 62, and 63 which are the 
main focus of interest here. Similarly, the range of N used 
in our fits was 100 < iV < 300. In Fig. H we show a2(i?), 
as^R), 62 (i?), and 63 (i?), fitted via Eqs. dSS]) and ([Ml), for 
several values of R spanning four orders of magnitude. 
The lines depict the anisotropy factors from Eq. ([Ml) , 
namely 52 = (2siii)~^ (dashed) and 53 = (2sit^)^-^ (fuU) 
multiplied by the pertinent values of 02(1) and 62(1) [for 
52 ] or 03(1) and 63(1) [for 33 ]. Once again, the agreement 
is perfect. The only case for which the higher-order terms 
{tti or 64) made any perceptible difference was for a^lR) 
at R> 30. 

Our results provide direct numerical evidence that the 
coefficient ratios 62/02 = 24.06(2) and 63/03 = 19.3(6) 



For fc = 1, this agrees with the conformal-invariance 
results of Eqs. © and ©; for fc 3, Eq. ^ gives 
63/03 = -240/7 [28[. 

Our results above for classical Ising spins differ from 
those for the QIC in that: (i) both even and odd powers of 
occur, in free-energy as well as in correlation-length 
expansions; and (ii) 63/03 = -1-19.3(6), incompatible with 
the value given in Ref. [2^ 

Note however that, when one considers the two- 
dimensional classical Ising model with the TM running 
along the diagonal Q, in the corresponding version of 
FBC, a picture closer to that found for PBC emerges, 
namely [23[: only odd powers of A^^^ occur in Eqs. ^ 
and (O, and the value 63/03 = —240/7 is reproduced. 



IV. S=l 

We considered 5* = 1 Ising systems on a square lattice, 
with both PBC and FBC, and isotropic couplings only. 
The critical temperature is known rather accurately [30| , 
J/keTc = Kc = 0.590473(5). 

In this case, no closed-form expressions for the TM 
eigenvalues are forthcoming, so one must rely on nu- 
merical diagonalization. The first consequence of this 
fact is that the assumption of only integer powers in 
Eqs. ^ and ^ must be reanalyzed. Indeed, while in 
Sections |ll] and IIIII one could verify directly from the re- 
spective closed-form equations that no noninteger powers 
of A^~^ were allowed, here this possibility does not arise. 
Furthermore, it has been shown for models very closely 
related to the standard Ising model that fractional pow- 
ers occur in corrections to scaling [31, 32]. It was con- 
jectured that these would take the form A^^/"^, clearly 
a very important term in the current context. However, 
for the S = 1 Ising model on a square lattice, it has been 
numerically shown that the amplitude of a hypothetical 
N~^/^ term is most likely zero [32], so in the current 
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Section at least, one can retain Eqs. ^ and ^ in their 
original form. 

Secondly, the range of strip widths within practical 
reach is much restricted in comparison with S = 1/2 
systems. We used 4 < iV < 16. Such a narrow range 
was, by far, the most quantitatively relevant source of 
systematic inaccuracies in our estimates of corrections to 
scaling, far outweighing, e.g., the uncertainties in Tc. 

In order to assess the associated effects, we produced 
fits of free-energy and correlation-length data for sets of 
S' = 1/2 data restricted to the same range of N. For 
PBC, we took truncated forms of Eqs. @ and ([S]) using 
the exact values of /oo, oi, and 6i, with {ak,bk} as ad- 
justable parameters for fc = 3, 5, 7, and zero otherwise. 
The k = 7 terms were included in order to increase sta- 
bility for the fc = 3 and 5 ones. By further restricting the 
range of data fitted to 10 < < 16, we found very good 
agreement with the known values of 03, &3, while 

for 05 and 65 deviations were of order 5% (see Table 

Turning to 5 = 1 with PBC, allowing for 02, ^2 7^ 
in Eqs. ([4]) and ([5]) gave fitted values of order 10"'^ — 
10~^ (compared with 03, 63 of order 10~^). We take 
this as signalling that, very likely, 02 = 62 = 0. Taking 
04, bi 7^ produced uncertainties of 50% or more in 
the corresponding estimates. This latter fact does not 
provide as compelling an argument to assume 04 = 64 = 
as the preceding one for 02, 62. However, in view of the 
limited number of data available for fitting, we decided 
that this was the most prudent route to take. 

Using 10 < iV < 16 and proceeding as described above 
for S — 1/2. we found the results shown in the last col- 
umn of Table ID Even assuming the systematic error in 
this case to be two orders of magnitude larger than that 
for S — 1/2, one gets 63/03 — 1.50(6), still at least 10 er- 
ror bars away from encompassing the S = 1/2 value. We 
refrain from attaching much significance to the estimates 
of b^/a^, due to the large uncertainty in 65. 

For S" = 1 systems with FBC, we fixed oi = 7r/48 Q, 
61 = 7t/2 [7]. Upon extrapolation of both PBC and 
FBC data, the non-universal bulk free energy is esti- 
mated as /oo = 1.317600(1). The surface free energy is 
/surf — —0.095187(1). Although the latter quantities are 
immediate byproducts of TM calculations, their value for 
S — 1 Ising spins on a square lattice does not seem to be 
available in the published literature Free-energy fits 
assuming 02, a3, and 04 as free parameters (the latter, for 
the purpose of stabilization of the former two), afc = for 
k > 5, gave a2 = -0.0238(2), 03 = 0.019(3), i.e., both 
of the same order of magnitude, contrary to the corre- 
sponding case for PBC. With similar assumptions for fits 
of correlation-length data, we obtained 62 = —0.5719(1), 
63 = -0.562(1). 



V. SECOND-NEIGHBOR COUPLINGS 

For square-lattice S = 1/2 spins with nearest- neighbor 
(next-nearest neighbor) couplings J (J'), we considered 



Table I. For 5 = 1/2 and 1, coefficients a^, bk from Eqs. @ 
and (O, and their ratios bk/ak- Columns 2 and 4: calculated 
from fits of free-energy and correlation-length data for strips 
with PBC across, and 10 < < 16, to truncated forms of 
those Equations (see text). Quoted uncertainties refer exclu- 
sively to the fitting procedures, i.e., no account is taken of 
likely systematic errors. Column 3: exact values from Ref. 0. 



Type 


S = 1/2, fit 


S ^1/2, exact 


S = 1, fit 


0.3 


0.15082(1) 


0.15073 . . . 


0.14772(2) 


as 


0.365(1) 


0.39214... 


0.313(3) 


63 


0.32305(1) 


0.322982 . . . 


0.2212(2) 


65 


0.766(3) 


0.79692 . . . 


0.11(5) 


bs/az 


2.1420(3) 


2.142857... 


1.497(2) 


brjas 


2.098(3) 


2.032258 . . . 


0.35(16) 



both interactions ferromagnetic and J' / J — 1. A gain , 
the critical point is known to excellent accuracy [3J|, 
= 0.1901926807(2). 

Once more, one must use numerical diagonalization of 
the TM since no closed-form expressions are available for 
the eigenvalues. We took 4 < A'' < 22, a significantly 
broader range than was feasible for S' = 1 in the preced- 
ing Section, but not in any way comparable to the leeway 
one has for S' = 1/2 with first- neighbor interactions only. 

Similarly to Section IIVI one must investigate whether 
noninteger powers show up in the corrections to scaling, 
N-^/^ being a likely candidate [Ulill. We did this by 
fitting our PBC free-energy and correlation length data 
respectively to: 



In — fa 







127V2 ^ 




TT 


bxe 







(38) 



where the adjustable powers Xf,x^ represent the domi- 
nant non-universal corrections. From fits of data in the 
range [iVo,22], we found Xf = 3.83(2), 3.971(3), 3.981(3) 
respectively for Nq — 4, 12, and 16, and x^ = 2.78(3), 
2.921(6), and 2.937(3) for the same sequence of Nq. So 
it is apparent that a;/ — > 4, a;^ — > 3 with increasing N. 
Comparing with Eqs. (|4]) and ([5]), we conclude for the 
absence of fractional powers such as iV"*/'^ here. 

For strips with PBC across, our analysis was then con- 
ducted along the lines described for = 1 in Section HVl 
Contrary to the S = 1 case, allowing for aj, bj did 
not improve stability of lower-order coefficients, and we 
decided to keep both to zero. The optimum range of 
widths for our fits was now 15 < < 22. We found 
03 = -0.09626(6), 05 = 0.210(4); 63 = -0.3305(3), 
65 — 1.96(1). From this, we estimate 63/03 — 3.43(1), 
which is again at variance with the S = 1/2 value Q 
15/7 = 2.14286.... 

For FBC, the known universal coefficients are oi = 
71-/48 y, bi = n/2 0. Combining PBC and FBC 
data, the extrapolated free energy per site is foe = 
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0.82926462(1), while the surface free energy is fsuri = 
—0.0895385(1). Estimates for these quantities are not 
quoted in pubhshed work on the next-nearest-neighbor 
S = 1/2 Ising model using TM techniques jH, [H. We 
attempted free-energy fits, at first using 02, a^, and 04 
as free parameters, and = for fc > 5. Similarly 
to the PBC case, allowing 04 to vary did not improve 
stabilization of 02 or 03, so we set 04 = 0. We thus 
found 02 = 0.00994(2), ag = -0.0096(1). From fits of 
correlation-length data, we obtained 62 = 0.23236(3), 
&3 = -0.529(1). 

VI. DISCUSSION AND CONCLUSIONS 

We have examined subdominant corrections to scaling 
for critical Ising systems on strip geometries. One of our 
main goals has been to check the extent to which the 
constant value of coefficient ratios, expressed in Eq. (|11L 
remains valid within the broader Ising universality class. 

In Section [n] we considered Ising S — 1/2 systems, 
on strips with PBC across. We investigated the effects 
of anisotropic interactions, extending the framework in- 
troduced in Refs. [l^ and [l^ and providing numerical 
evidence that the non-universal coefficients 03 and 63 of 
Eqs. Q and ^ indeed follow the prediction given by 
Eq. ((T5|) . As a byproduct, the validity of Eq. (fTTj) has 
been directly verified within four orders of magnitude of 
anisotropy variation for this case. 

In Section Hill for strips of spin-1/2 systems with EEC 
along one of the coordinate axes, we examined ways in 
which the non-constant frequency spacing in the eigen- 
value spectrum can be dealt with, in order to make the 
sum in Eq. (EU)) amenable to treatment via the Euler- 
Maclaurin summation formula. The lowest-order terms 
of the resulting expansion are shown to agree with known 
results. 

From the correlation- length expression, Eq. (HI]), we 
showed directly that both odd and even powers of inverse 
strip width are expected in corrections to scaling, and ex- 
plicitly evaluated the two lowest-order non-universal co- 
efficients [ see Eq. ((55)) ] . Generalization to anisotropic 
systems is given in Eq. ([M)) . where one can see that 
the first- and third order anisotropy factors (respectively, 
1/si and l/2sitl) are the same as those for PBC [see 
Eqs. (ini), (dSl), and 

We also found numerically that the amplitude ratios 
bk/cik remain constant, for k = 2 and 3, upon introduc- 
tion of anisotropic couplings. 

Sections ll VI and IVl deal respectively with S — 1 systems 
with first-neighbor interactions, and spin-1/2 ones with 



both first- and second-neighbor couplings. For PBC we 
find that, in both cases, the ratio 63/03 differs consider- 
ably from the value 15/7 = 2.142857 . . . found in Ref. H 
for S — 1/2, first-neighbor couplings only. We quote 
63/03 — 1.50(6) for the former, and 3.43(1) for the latter. 
For FBC, comparison of bk/o-k ratios with those pertain- 
ing to S = 1/2 systems gives 62/02 = 24.0(2) {S = 1), 
23.4(1) (next-nearest-neighbor), 24.06(2) {S = 1/2, first- 
neighbor). Although the error bars do not quite overlap, 
it appears that a constant value of this ratio cannot be 
definitely discarded. However, no such regularity is seen 
for 63/03, its value being respectively —30(5), 4-55(1), 
and -1-19.3(6) in each case. 

Overall, it seems that both even and odd powers of 
always show up in Eqs. (|4]) and ([5]), for critical 
Ising strips with FBC along one coordinate axis. On the 
other hand, for PBC only odd ones occur. Concurring 
remarks can be found in the literature [26|: however, 
it seems difficult to prove such a statement rigorously. 
So far, one has to rely on case- by-case analyses, as was 
done here. As pointed out at the end of Section lllli 
considering the version of FBC with the TM running 
along the diagonal [3| is enough to restore a picture very 
similar to that holding for PBC Thus, the behavior 
of subdominant corrections to scaling is sensitive to what 
might appear to be a minor technical detail. 

The constant value of amplitude ratios is maintained 
upon varying anisotropy for 5* = 1/2 systems with first- 
neighbor couplings, either with PBC or FBC; however, 
it does not seem to survive changes in spin 5, or intro- 
duction of further neighbor interactions. We have thus 
established that the observed, apparently universal, con- 
stant amplitude ratios pertain to a limited subset of sys- 
tems which are in the broader Ising universality class. 
It remains to be further investigated whether the close 
values found for 62 /02 with FBC in the three cases are 
indeed an indication of an actual constant ratio. 
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